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Abstract. We present sub-arcsecond, Multi-Element Radio Linked Interferometer (MERLIN) observations of the decimetre 
radio continuum structure and neutral hydrogen (Hi) absorption from the nuclear region of the starburst galaxy NGC4194 (the 
Medusa Merger). The continuum structure of the central kiloparsec of the Medusa has been imaged, revealing a pair of compact 
radio components surrounded by more diffuse, weak radio emission. Using the constraints provided by these observations and 
those within the literature we conclude that the majority of this radio emission is related to the ongoing star-formation in this 
merger system. 

With these observations we also trace deep Hi absorption across the detected radio continuum structure. The absorbing 
Hi gas structure exhibits large variations in column densities. The largest column densities are found toward the south of the 
nuclear radio continuum, co-spatial with both a nuclear dust lane and peaks in 12 CO(l— >0) emission. The dynamics of the 
Hi absorption, which are consistent with lower resolution 12 CO emission observations, trace a shallow north-south velocity 
gradient of ~320 km s" 1 kpc~' . This gradient is interpreted as part of a rotating gas structure within the nuclear region. The Hi 
and CO velocity structure, in conjunction with the observed gas column densities and distribution, is further discussed in the 
context of the fuelling and gas physics of the ongoing starburst within the centre of this merger. 

Key words, galaxies: individual NGC4194: hydrogen absorption- galaxies: Radio Lines -galaxies: active - Star-formation: 
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1. Introduction 

The Medusa merger (NGC4194) is a nearby (D=39Mpc) star- 
forming galaxy of intermediate infrared luminosity (Lfir = 
8.5 x 10 10 L Q ). As such the Medusa merger falls into a lumi- 
nosity class, considerably lower than the well-studied ultra- 
luminous infrared galaxies (ULIRGs) but it still appears to 
show areas of intense starburst activity (Armus, Heckman & 
Miley 1990; Prestwich, Joseph & Wright 1994). The optical 
morphology of the Medusa merger, as its name suggests, traces 
clumpy and extended features which stretch ~ 1' to the north of 
the main optical nucleus. These 'hair-like' extensions, indica- 
tive of tidal tail features resultant from the galaxy's ongoing 
merging event, contain significant volumes of molecular gas 
(Aalto, Huttemeister & Polatidis 2001) although no indications 
of ongoing star formation in the tidal tail have, as yet, been 
detected (Armus et al. 1990). The majority of star formation 
in the Medusa merger is concentrated in the extended nuclear 
starburst (Armus et al. 1990). 

A detailed study of the molecular gas (CO J=l— >0) com- 
position and kinematics in NGC4194 has been previously 
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made at arcsec angular resolutions at the Owens Valley Radio 
Observatory (OVRO) by Aalto & Huttemeister (2000, hereafter 
AH00). This investigation revealed that the 12 CO emission is 
extended over an area of 25"(4.7kpc) covering the central and 
north-eastern part of the optical galaxy. The CO emission also 
traces the two prominent dust lanes that cross the central region 
and extends into the northern tidal tail. The majority of the CO 
(~70%) is found with the central 2 kpc of the galaxy, with 15% 
of this gas residing in a compact region 1'.'5 south of the radio 
nucleus. 

We have used the UK's Multi-Element Radio Linked 
Interferometric Network (MERLIN) to observe both the 
1 .4 GHz radio continuum and neutral hydrogen (Hi) absorption 
structure of the central few hundred parsecs of NGC4194 at 
sub-arcsec angular resolutions. This constitutes the highest an- 
gular resolution study made to date of both the radio continuum 
and neutral gas structure in this merger galaxy. In this paper we 
present the results of these observations and place them in con- 
text with other observations of the radio continuum, and neu- 
tral and molecular gas in NGC 4194. The primary goal of this 
work is to provide high linear resolution information about the 
neutral gas within the core region of this star-forming galaxy, 
enabling the composition and dynamics of the central region 



2 



Beswick et al.: Hi absorption in the Medusa Merger 



to be investigated. In conjunction with lower resolution data 
(e.g. AHOO; Aalto et al. 2001) these observations can be used 
to more fully chart the gas dynamics and fuelling of this inter- 
mediate luminosity starburst galaxy. 

This paper is organised as follows. The observations and 
the data reduction processes undertaken are described in detail 
in the next section. This is followed by the presentation of the 
radio continuum and Hi absorption results. In the subsequent 
sections we will discuss, in detail, these observational results, 
initially drawing conclusions regarding the neutral hydrogen 
and radio continuum structure of the central 0.5 kiloparsec re- 
gion of the Medusa and then by placing these results in the 
context of other observations of this source's gas content. 

Throughout this paper we will assume a distance of 39 Mpc 
to NGC4194 (Hf^Skms-'Mpc- 1 ), implying that at the dis- 
tance of the source 189pc subtends an angle of 1 second of 
arc. 

2. Observations and Data Reduction 

2.1. MERLIN 

NGC4194 was observed using the MERLIN array, includ- 
ing the 76-m Lovell Telescope (Thomasson 1986), on 2003 
March 15 at 1408.66 MHz. NGC4194 was observed for 
~18h interspersed with regular observations of the nearby 
phase calibration source 1205+545. Throughout the observa- 
tion dual hands of circular polarisation were recorded over 
an 8-MHz band which were correlated into 64 spectral chan- 
nels, each of bandwidth 125 kHz providing a velocity resolu- 
tion of ~26kms~ 1 channel -1 , centred upon the rest velocity of 
NGC4194, 2560kms _1 . A single observation of OQ208 was 
used to obtain passband solutions along with a 30min inte- 
gration on 3C 286, which was assumed to have a flux density 
of 14.794 Jy at 1408 MHz (Baars et al. 1977) and was subse- 
quently used to calibrate the flux density scale for all sources. 

Initial editing and calibration of these NGC4194 data 
were carried out using local MERLIN software at Jodrell 
Bank Observatory, prior to these data being read into the 
Astronomical Image Processing Software (aips). Within aips 
further editing and bandpass calibration were applied to the 
data set. Phase corrections derived from the phase calibrator 
source 1205+545 were also applied to these data, along with 
further phase corrections derived from self-calibration upon 
the continuum of NGC4194. The calibrated uv data set was 
Fourier-transformed to produce spectral line data cubes with 
various tapering and data weighting schemes used in order to 
produce images that were sensitive to diffuse large scale radio 
emission and high resolution images of the compact compo- 
nents. The line-free channels were combined to produce con- 
tinuum images which were subsequently subtracted from the 
spectral-line cubes. These continuum images along with the 
continuum subtracted spectral-line cubes were deconvolved us- 
ing a clean based algorithm to remove the instrumental re- 
sponse (Hogbom 1974). The cleaned continuum images and 
continuum-subtracted cubes were then recombined and further 
analysis of these data were undertaken using standard spectral 
line tasks within aips. 
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Fig. 2. Uniform weighted MERLIN 1.4 GHz radio contin- 
uum image of the core of the Medusa merger with a synthe- 
sised restoring beam of 0717 x 0/15. Contour levels start at 
0.35 mJy bm~' and increase by factors of V2. 



These MERLIN data have a shortest projected baseline 
of ~13km between the Tab ley and Jodrell Bank telescopes 
and hence these data are not sensitive to extended compo- 
nents greater than a few arcsecs. The consequences of this 
data-set including few low-order spacings is a hole in the uv 
coverage toward the centre of the uv plane. This feature of 
sparsely filled interferometric arrays, such as MERLIN, causes 
an under-sampling of the low-order spacings and hence an in- 
sensitivity to regions of smooth, extended flux. This effect can 
result in a bowl-shaped reduction in the flux density toward 
extended sources (see, for example Wills, Pedlar & Muxlow 
1998). This not only reduces the flux density of radio con- 
tinuum components derived from images but also causes the 
calculated Hi absorption opacities to be over estimated due to 
the reduction in the value of the extended continuum emission. 
The total flux measured from these observations is 28 mJy com- 
pared to a lower resolution (beam size of 45") VLA 1.4 GHz 
flux of 91 mJy (Condon et al. 1990). The distribution of the 
radio emission detected in this VLA image covers an area is 
approximately 4 times larger than the region shown in Fig. 1 . 
This implies that these observations are spatially filtering out 
approximately two thirds of the 1 .4 GHz flux. If it is assumed 
that this missing flux is spread over a Gaussian with an area 
matching the component detected by Condon et al., the maxi- 
mum correction to the observed MERLIN map that is derived 
is 0.05 mJy at any single point in the image. In the case of the 
observations presented in this paper we have further tested for 
the significance of these missing low-order spacings by taking 
slices across the images in several directions. Although some 
evidence for a slight reduction in the flux density of the im- 
ages produced was found, it has not been directly corrected for 
in the images and spectra present since the correction required 
was found to be significantly smaller than the channel based 
noise levels in these data. 
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Fig. 1. Overlay of uv-tapered (800 k/l) 1.4 GHz radio continuum emission in green contours and OVRO 12 CO (AH00) in pale 
blue contours over a F606W WFPC2 HST image. The contour levels on the radio continuum image are -1, 1, 2, 4, 8, 16 and 32 
times 0.3 mJy birr 1 . The '"CO contours are 1, 2, 4, 8, 16 and 32 times 0.5 Jy beam -1 kms -1 . The false-colour scale of the optical 
image is linear and arbitrary. The radio and CO images have synthesised beams of 07446x0'.'320 and 2'.'06xr.'70 respectively. 
The astrometrical alignment of the HST image with the other two images has an rms error of ~ 1 00 mas and is described in section 
2.2. Labels a, b, c, d and e indicate the 12 CO peaks as observed by AH00. 



2.2. Hubble Space Telescope (HST) 

In order to afford a good comparison between the sub-arcsec 
radio observations presented in this paper and the optical emis- 
sion from this galaxy archival HST observations have been 
acquired. These NGC4194 data (U2E66Y01T) consisted of 
single 600 s integration using a broad-band F606W filter (see 
Fig. 1). Initially these data were reduced using both stsdas 
commands within iraf and starlink software before being ex- 
ported as a fits file. The procedure of applying accurate astrom- 
etry to the HST field was achieved by assigning known coordi- 
nates, obtained from the US Naval Observatory star catalogue, 
to a number of stars spread throughout the WFPC-2 field. This 
procedure resulted in a linear transformation and rotation be- 



ing applied to the HST image in-order to minimise the residual 
position errors of the known star positions. No shear was ap- 
plied to the image during this process. It is estimated that the 
absolute positional accuracy of the aligned MERLIN and HST 
images is <0.1 arcsec. 

3. Results 

3.1. The 1 .4 GHz radio continuum 

A uv-tapered (800 k/l) 1.4 GHz image of the central region 
of NGC4194, derived from the line-free channels, is shown 
as contours in Fig. 1, overlayed upon a false-colour archival 
WFPC2 HST image. A u v-taper was applied to these data in or- 
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Table 1. The 1.4 GHz radio continuum properties of the core components. 3 
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a Listed properties are for the two compact nuclear components. Flux densities and positions have been derived from Gaussian fits to our 
uniformly weighted image (Fig 2). 



der to produce images which are as sensitive as possible to the 
weak and diffuse extended radio emission. A measured noise 
level of 0.1 mJybrrr 1 has been determined for Fig. 1 and the 
image has an angular resolution of 0745x0732. Figure 2 shows 
a high angular resolution, uniformly weighted, contour image 
of the core region of the galaxy. This figure has an angular 
resolution of 07 17x0'.' 15 and a noise level of ~0.1 mJybrrT 1 . 
Table 1 lists the 1 .4 GHz flux densities and Gaussian fitted sizes 
of the two, bright, compact components seen in the uniformly 
weighted image (Fig. 2). The radio structure of the central few 
hundred parsecs of NGC4194 is dominated by a pair of com- 
pact radio components (Fig. 2). These two components are sep- 
arated by -0735 (-65 pc) along a PA of 120±5°. In the uv- 
tapered image shown in Fig. 1, and enlarged as the central im- 
age of Fig. 3, it can be seen that these two compact radio com- 
ponents are only partially resolved and are embedded in a halo 
of diffuse radio continuum which extends in a southerly direc- 
tion. 

3.2. Neutral hydrogen absorption 

As can be seen in Fig. 3, extensive Hi absorption is detected 
against the majority of the background radio continuum in 
the uv-tapered MERLIN image. Similarly Hi absorption is 
also found against the core region in spectral line cubes made 
with uniform weighting. The Hi absorption properties derived 
from the uv-tapered spectral line cube are listed in Table 2. 
These properties include the central positions of the spectra 
shown in Fig. 3, in addition to the Gaussian fitted velocities 
and widths of the absorption lines. Table 2 also lists derived 
properties of these absorption lines, such as their peak opac- 
ities and column densities of Hi along these lines of sight. 
The column densities listed have been calculated assuming that 
N H = 1.823 x 10 18 T spin /rdV. Although T spin can be reason- 
ably assumed to be 100K (Maloney, Hollenbach & Tielens 
1996) it should be noted that the value of T sp ; n is dependent 
upon the physical conditions of the gas and will therefore vary 
significantly especially in the extreme environment of a nuclear 
starburst. These considerations will be discussed in more detail 
in a latter section. 

In Fig. 4 we present three grey-scale images of the centre of 
the Medusa presenting the Hi absorption detected against the 
radio continuum. The left-most (Fig. 4a) of these panels shows 
a grey-scale map of the Hi absorption toward areas with signifi- 



cant radio continuum background. As can be seen in this image 
the line signal varies, with the absorption deepest toward the 
eastern edge of the brightest radio component. The map of in- 
tegrated absorption shown in Fig. 4a does not represent the true 
Hi column density since it is dependent upon the background 
continuum strength. In the central panel (Fig. 4b) the distribu- 
tion of the integrated Hi opacity is shown. This panel shows the 
distribution of neutral hydrogen and is directly proportional to 
the Hi column density. It should, however, be noted that the pro- 
cess of deriving maps of optical depth involves the division of 
spectral line data-cubes with their associated radio continuum 
images, resulting in the opacity images produced being suscep- 
tible to errors in areas where the radio continuum is small (i.e. 
toward the edges of the image). Measures to minimise this ef- 
fect have been made, although some unphysical values in a few 
pixels toward the edges of the image may remain. However 
considering this factor the variations in opacity seen are signif- 
icantly larger than the errors in the image. Plus the values and 
features seen in this image are consistent with those observed 
and calculated from the absorption spectra (see Fig. 3) and ab- 
sorption maps (Fig. 4a). The Hi integrated opacity shows an 
increase in the column densities south of the radio continuum 
peak which is consistent with the positions of the peaks of the 
CO emission observed by AH00 (labelled in Fig. 1 & 4a). 

Figure 4c shows that the velocity structure of the Hi fol- 
lows a distinct north to south velocity gradient on the largest 
scales along with an N-W to S-E velocity gradient across the 
two bright, compact radio components. Fig. 5 shows a posi- 
tion (Declination) verses velocity contour plot of the absorption 
averaged over the entire R. A. extent of the radio continuum 
source. This figure complements the two dimensional veloc- 
ity distribution shown in Fig. 4c by showing in detail the north 
to south velocity gradient in the source. As can be seen from 
these figures the velocity structure illustrated in Fig. 4c traces a 
gradient of ~60 km s~ 1 arcsec~ 1 . 

4. Discussion 

4.1. High resolution radio continuum structure of the 
centre of the Medusa 

In the highest angular resolution (07 17x0715) image produced 
using these MERLIN data (Fig. 2) the majority of the radio 
structure of the Medusa galaxy is resolved away leaving a 
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Fig. 3. Montage of 800-k/l uv-tapered radio continuum structure and Hi absorption spectra. Contour levels of the central image 
start at 0.3 mJybm -1 and increase by factors of V2 with an angular resolution of 0745 x 0732. 

pair of compact nuclear components (see Table 1 & Fig. 2) the core regions of many low powered Seyfert galaxies (e.g. 
with a small, weak extension toward the north-east of the Thean et al. 2000). However, the lower angular resolution uv- 
brightest component. This double structure is reminiscent of tapered images show these components to be sitting in an ex- 
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Fig. 4. a) The (leftmost) image shows the integrated Hi absorption signal across the centre of the Medusa merger. This has 
been derived from the continuum subtracted spectral line cube. The grey-scale flux density range is 0— > -185 mJy birr 1 km s _1 . 
Note the absorption shown is dependent on the strength of the continuum emission. The two symbols (a and b) shown on the 
left-most image represent two of the CO emission peaks from AHOO. A single contour 1 mJybirT 1 radio continuum contour is 
overlayed. Areas external to this contour have the lowest levels of radio continuum and have been excluded in the derivation 
of Fig 4b. b) The (central) image shows the integrated optical depth ( JrdV) across the same region as a. The grey-scale range 
shown is 15— >500 and the contour levels are 10, 20, 40, 80, 160, 320 and 640 times 100. It should be noted that errors in this 
image increase in areas where the radio continuum flux density is at a minimum. This effect results in some erroneous values 
toward the edges of the image, in order to remove the most severe of these effects this image has been blanked in areas where the 
radio continuum is less than 1 mJybeam -1 . c) The (right-hand) map represents the Hi absorption velocity field over areas where 
significant radio continuum and absorption line signal were detected. The grey scale range of the image is 2430— > 26 10 km s -1 , 
the overlayed contours have increments of 25 km s _1 between 2475-> 2650km s _1 . 



tensive, diffuse, radio emission toward the south of these com- 
pact components (see Figs. 1 and 3). These uv-tapered im- 
ages (Fig. 1 and 3) are consistent with other studies of this 
source previously made at higher frequencies using the Very 
Large Array (VLA) (e.g. Ulvestad, Wilson & Sramek 1981). 
Fig. 3 is still insensitive to many of the weak, larger scale fea- 
tures observed by Ulvestad et al. (1981). More recent multi- 
frequency (5 to 15 GHz) VLA observations (S. Neff, private 
comm.; Neff, Campion & Ulvestad 2001) show that in addi- 
tion to the compact components and diffuse emission observed 
in these MERLIN 1.4 GHz data more extended emission is 
present. In particular, in deep 5 GHz VLA images (S. Neff, 
private comm.), an 'arm-like' extension arcing from the north 
of the nuclear region and ending approximately 4 arcsec to the 
west of this region in a strong compact component, is observed. 
This feature is also seen in VLA images published by Ulvestad 
et al. (1981). Although some traces of this 'arm-like' emission 
can be seen in our uv-tapered images (see Fig. 1) the majority 
of these features are not detected by our observations. Neff et 
al. (2001) also report -50 compact sources (> 5x 10 18 W Hz" 1 ). 
Of the brightest of these compact sources ~50% nominally 
have flat spectral indices and maybe associated with thermal 
radio emission from Hn regions, whereas the remainder have 



steeper spectral indices that are indicative of emission dom- 
inated by non-thermal mechanism (such as supernovae rem- 
nants or AGN). Our sub-arcsec 1.4 GHz observations are, as 
mentioned above, insensitive to much of the diffuse extended 
radio emission detected with the VLA, but more sensitive to 
emission from synchrotron sources than thermal sources; the 
combination of these effects and our sensitivity levels accounts 
for the differences in observed emission structures. 

4.2. The neutral hydrogen absorption 

High optical depth (Tp^ >0.5) neutral hydrogen absorption 
is detected against almost all of the radio continuum emission 
detected (Fig. 1). The absorbing Hi gas is observed to show 
significant variations in opacity across the radio continuum 
source. The Hi column densities observed against the radio 
continuum emission are consistently higher toward the south- 
ern half of the nuclear radio continuum. These variations can 
be seen in both the maps of the integrated Hi absorption and 
opacity (Fig. 4a, b) as well as in the individual spectra (Fig. 3 
and Table 2). This implies that the Hi observed, is not part of 
a uniform homogeneous screen of neutral foreground gas. It is 
likely that the large scale variations in the Hi opacities observed 
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are the result of high gas column densities along the lines 
of sight of other large scale features, such as the dust lanes. 
Whereas smaller, more localised, variations seen in Fig. 4b (on 
size scales of order a beam [~50-60pc]) are likely to be due 
to single, or the superposition of several, gas clouds along the 
sight line. In Fig. 1 the distribution of the radio continuum rel- 
ative to the WFPC2 R-band image shows a thick dust lane arc 
that crosses the southern part of the MERLIN radio continuum 
in an approximately east-west direction. The position of this 
dust lane, which is also co-spatial with the CO emission (red 
contours in Fig. 1; AHOO), is coincident with the area of high- 
est opacity Hi absorption. It is also reasonable to assume that 
the gas, both CO and Hi, is associated with this dust lane (the 
relative distributions of these constituents is discussed in more 
detail in Sec. 4.4.1). Hence this implies that the radio contin- 
uum originates from deeper within the source than these com- 
ponents or that these elements of the neutral ISM and the radio 
emitting material are mixed resulting in these observations un- 
derestimating the true Hi column densities. 

The Hi velocity field shown in Fig. 4c traces a dominant 
but shallow velocity gradient with a north to south direction. 
This gradient is also shown in Fig. 5. In this averaged position- 
velocity slice the broad absorption tracks a velocity gradient of 
~60kms" 1 arcsec^ 1 (~320kms~ 1 kpc~ 1 ). This velocity gradient 
extends over ~2'.'5 (470 pc) of the background radio continuum 
source. The north-south orientation of the velocity field is con- 
sistent with high resolution CO observations from the same re- 
gion (AHOO). The details of the velocity field (Fig. 4c) show 
that the primary gradient deviates from a north to south direc- 
tion against the brighter northern radio continuum components. 
Hints of this directional change in the Hi velocity gradient are 
visible in the highest angular resolution (~T.'7) CO data pre- 
sented by AHOO. 



4.3. The source of the radio emission in the Medusa 

Although the Medusa merger is an order of magnitude fainter 
at infrared wavelengths than the well studied ULIRG galaxies, 
such as Arp 220, it has in the recent past undergone, and contin- 
ues to undergo, significant star-formation activity. This activity 
has presumably been triggered by its advanced merger state. 
However, it is also widely accepted that gas rich mergers, such 
as NGC 4194, can also host AGN activity in the remnant nuclei 
of the merged sources (Genzel et al. 1998). This phenomenon 
is observed in some composite starburst-AGN sources (e.g. 
Beswick et al. 2001; Gallimore & Beswick 2004; Kewley et al. 
2000) and is also evident in some minor mergers which host 
more powerful AGN such as radio galaxies (Heckman et al 
1986). As discussed previously the radio continuum observed 
from the Medusa merger, at 1 .4 GHz, traces both the areas of 
starburst activity, as well the more compact radio components 
that may be AGN related. In the rest of this section we will 
discuss both the nuclear starburst of the Medusa and its rate of 
star-formation, as well as briefly explore the possible evidence 
that a weak AGN exist in the core of this galaxy. 



4.3.1 . An AGN component? 

NGC 4194 is optically classified as a peculiar, blue compact 
galaxy with a strong Hn emission line spectra and a blue contin- 
uum indicative of a young stellar population (Liu & Kennicutt 
1995). However, as discussed above the two compact radio 
components imaged in the centre of this galaxy show a struc- 
ture that is similar to that observed in many low luminosity 
Seyfert galaxies. These two components are relatively low in 
luminosity, 5.48xl0 20 and 2.22x1 20 WHz but are compact 
on MERLIN scales (<25pc), implying that the strongest of 
these components has a brightness temperature in excess of 
4.5xl0 4 K. Searches for an AGN-like radio core using VLBI 
techniques in the Medusa have so far not detected any bright 
compact radio components (Lonsdale, Lonsdale & Smith 1992; 
vanBreugel et al. 1981). 

In addition to this continuum information these spectral line 
observations, along with those of AHOO provide information 
regarding the gas velocity distribution within the central region 
of NGC 4 194 which can be used to place some limits on the 
mass enclosed within the central few arcsec. From these Hi ob- 
servations we observe a velocity gradient in the neutral gas of 
~320kms~ 1 kpc~ 1 over the central 275 (470pc) which, if as- 
sumed to be in circular motion, implies a limit on the enclosed 
mass of ~10 6 M . Such a limit constrains the mass of any black 
hole within the centre of the NGC 4194. 

However, considering the limits on the luminosities and 
sizes of these radio continuum components, the presence or 
not, of a weak AGN located at the position of one of these 
two components cannot be categorically ascertained. The ev- 
idence from the VLBI studies of Lonsdale et al. (1992) show 
that no high brightness temperature synchrotron >1.8mJy at 
A =18 cm is detected. This supports the case for little or no 
AGN related radio emission. Although it should be noted that 
southern nucleus of the binary AGN NGC 6240 contains a very 
weak compact VLBA core with flux of »1 mJy (Gallimore & 
Beswick 2004) which is free-free absorbed at 18 cm. Further 
high resolution multi-frequency radio imaging and X-ray ob- 
servations are required to investigate the possible scenario that 
there is an AGN buried in the centre of this source. 

4.3.2. The starburst & star-formation rate 

As has already been stated the Medusa merger is dominated, 
in terms of its optical and infrared emission by starburst ac- 
tivity. Starburst emission results in relativistic charged parti- 
cles accelerating in supernovae remnants (SNRs) and radio su- 
pernovae (RSNe), which produces the radio emission, such as 
found in other nearby starburst galaxies (e.g. M82 - Muxlow et 
al. 1994). 

It has been shown, for example by Condon (1992) and 
Cram et al. (1998), that the 1.4GHz radio emission from star 
forming galaxies can be considered as a direct tracer for the su- 
pernovae rate. If this is assumed, the star-formation rate (SFR) 
of the source can also be estimated, albeit with an initial mass 
function (IMF) adjusted to produce massive enough stars to 
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Table 2. The Hi absorption properties of the central region of the Medusa merger. 
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AV(kms-') {&) 




Nf, xl0 19 T spin 
(atoms cm -2 ) 


l b 


402 


16V085 










2 b 


401 


16^545 










3 


4P.22 


16V206 


2482 


41 


0.115 


0.08 


4 


41 s .25 


16'.' 059 


2520 


49 


0.533 


2.02 


5 


41 s .23 


15'.'635 


2528 


58 


0.964 


11.27 


6 


4P.22 


14'.'915 


2587 


42 


1.350 


4.82 


7 


41*25 


14'.' 600 


2586 


52 


1.425 


5.89 



a Labels and positions listed refer to the centre of the areas over which the spectra in Fig. 3 are averaged. 

b No limits upon the optical depth are given in areas where the radio continuum flux density is too low to provide a physically important 
limit. 

Column densities have been calculated assuming Nh = 1.823 x 10 18 T sp ; n f TdV. These are derived from the absorption opacities and line 
widths over areas of approximately a beam in size. The value of T spin is assumed for the purposes of the discussion to be 100 K, however 
this variable will be discussed in more detail in Sec. 4.4.3. 



form supernovae that result in significant radio emission (i.e. 
M>8 M ) such that, 



SFR 



1.4 



5.5 x 



-1.4 



4.6 x 10 21 WHz 



— M yr 



-l 



(1) 



(e.g. Haarsma et al. 2000; Hopkins, Schulte-Ladbeck & 
Drozdovsky 2002). In this equation the multiplying factor of 
5.5 is included to convert the mass range 5-100 M used by 
Condon (1992) to the range 0. 1-100 M , assuming a Salpeter 
IMF. The low angular resolution (arcsecond), 1.4 GHz VLA 
flux of NGC4194 is 92.7 mJy (Condon et al. 1990), implies a 
total SFR1.4 «2OM yr _1 . In the case of the Medusa, and most 
starburst and ULIRG systems, the accuracy of the SFR as de- 
rived from optical data (e.g. Ha observations) will be severely 
affected by the large and variable nature of the obscuration in 
this source. Using Ha observations Storchi-Bergmann, Calzetti 
& Kinney (1994) quote a SFR Hq , «38M yr _1 (assuming 
H =50kms _1 Mpc -1 ) which corresponds to a more modest 
value of ~17M Q yr _1 , if H =75kms _1 Mpc _1 is assumed. This 
value shows a remarkable consistency with the SFRi 4 value 
obtained from the total radio flux considering the uncertainties 
that are introduced by the high levels of nuclear obscuration 
within NGC4194. 

In our higher resolution radio data only the emission from 
the intense, compact nuclear star-forming regions can be de- 
tected. Applying equation 1 to our total MERLIN flux of 
28mJy implies a SFR of ~6M yr 1 for the Medusa; com- 
pared to the far-infrared derived SFR/,>=6-7M yr _1 (AH00). 
Clearly in the case of this radio determination the value can 
only be considered as a lower limit for the total SFR, since it 
only accounts for star-formation in the compact nuclear region. 

Compact clumps of star-formation have been imaged in 
NGC4194 by Weistrop et al. (2004) at UV wavelengths using 
the HST. They find that there are ~40 knots within the centre 
of NGC4194. From their analysis they conclude that the total 
SFR within the youngest (<20Myr) of these knots, which are 
primarily located in the nuclear region in areas broadly consis- 



tent with the distribution of the 1 .4 GHz MERLIN radio con- 
tinuum, is ~6M yr _1 . This is consistent with our radio and 
far-infrared determinations of the SFR in the nuclear region. 

The structure and luminosity of the radio emission in the 
lower resolution 1 .4 GHz radio images (Fig. 3) can easily be re- 
produced via pure starburst mechanisms with only a relatively 
modest nuclear SFR. In this scenario there is no requirement 
for a compact radio AGN although it is not precluded. 

4.4. The ISM in the nucleus of the Medusa 

AAA . Distribution of Hi, CO and dust 

Using these absorption data only gas in front of the radio con- 
tinuum can be detected. However it is also apparent the Hi dis- 
tribution throughout the nuclear region is not an even screen 
of gas. There are significant variation in the measured opac- 
ity of the absorbing Hi (assuming a constant T^n 1 ), with high 
columns of gas situated toward the southern part of the radio 
emission. This area of increased Hi column density is spatially 
coincident with the increased optical extinction from dust that 
arises from the cross-nuclear dust lane (see Fig. 1). The impli- 
cation of this is that, in addition to the optical emission, the ra- 
dio emission is situated behind the dust lane and that it is prob- 
able that the obscuring dust and neutral gas are physically asso- 
ciated. This is commonly seen in many nearby active galaxies 
with nuclear dust lanes {e.g. Cole et al. 1998; Beswick et al. 
2002, 2003, 2004; Jackson et al. 2003). 

If we assume a standard Galactic gas-to-dust ratio (Stavely- 
Smith & Davies 1987), an inferred optical extinction can be 
derived from the Hi column density. Assuming these ratios the 
optical extinction at B-band is Ab = 8.62 X 10~ 22 NHmag. 
Consequently a column density of ~5xl0 21 atoms crrT 2 
implies an extinction along the same line of sight as the dust 



1 Variations in the value of T spin across the radio continuum will 



also effect the line strength (see Sec. 4.4.2). 
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lane (positions 5, 6 and 7), of A B «4-5 mag, assuming con- 
stant T S pi n »100K. This value is approximately 50 times larger 
than the implied levels of extinction against the compact radio 
components toward the north of this radio source which are not 
crossed by the nuclear dust lane (see Fig. 1). 

AH00 mapped the distribution of CO (J=l->0) in the 
nucleus of NGC4194 at high angular resolution using 
OVRO. In this study the CO emission in the velocity range 
2506-^2570 km s- 1 (see figure 7 of AH00) is also found to 
trace the primary dust lane of the Medusa. AH00 observed CO 
emission to have linewidths of FWHM ~150kms _1 , larger but 
still comparable to those observed here (see Table. 2). This is 
quite remarkable considering the angular resolution of these Hi 
observations is an order of magnitude higher than the CO emis- 
sion measurements. The consistency of these two linewidth 
measurements could imply that these two gas components, sit- 
uated along the same line of sight, are well mixed over a range 
of angular scales (~2— >0.2 arcsec). 

At the areas of strongest CO emission (a and b in Fig. 1 
and AH00) the CO gas surface density is -500-1000 Mope 2 
(~4-8xl0 21 molecules cm -2 ). These columns compare to the 
Hi absorbing column densities of ~6-12xl0 21 atoms cm -2 , as- 
suming T spin = 100K. These Hi columns are merely a lower 
limit since they are only sensitive to the gas in front of the ra- 
dio continuum. As stated above these observations also show 
an increased column density of Hi toward the radio continuum 
situated along the line of sight of the dust lane, which is also co- 
incident with the CO, in particular CO emission peaks a (AH00 
and Figs. 1 & 4). Thus it implies that the nuclear CO, Hi and 
dust are all situated along the same sight line. 

4.4.2. Dynamics of the cold ISM: the fuel for the 
evolving starburst 

Figures 4c and 5 show the velocity field and a north-south 
position-velocity diagram derived from these Hi absorption ob- 
servations. Comparing the sub-arcsec angular resolution Hi ab- 
sorption velocity field (Fig. 4c) with the centre of the ~1'.'7 an- 
gular resolution CO velocity field presented in Fig. 8b of AH00 
it is apparent that the 'global' nuclear CO and Hi dynamics are 
consistent. In both of these cases the general trend is for the 
gas to be rotating following a north-south velocity gradient of 
a few tens of kms 1 per arcsec; further confirming the likely 
association of these two components. 

On the scale of the OVRO CO emission the radio contin- 
uum peak is coincident with the kinematic centre of the galaxy. 
Centred close to this kinematic centre, these MERLIN obser- 
vations provide a 'close-up' view of the nuclear gas dynamics. 
In Fig. 4c the velocity field, although consistent with the CO, 
shows significant deviations and smaller scale velocity struc- 
ture that is smoothed away by the larger CO synthesised beam. 
In particular the isovelocity contours of Fig. 4c show a signifi- 
cant twist toward the west as they cross the area of peak radio 
flux density (see also spectra 3 and 4 of Fig. 3). This twist in 
the velocity field is hinted at in the CO observations of AH00. 
However AH00 conclude that the gas in this region undergoes 
approximately solid-body rotation. This is broadly confirmed 



2380 




-2.0 -1.5 -1.0 -0.5 0.0 

Offset position (arcsec) 

Fig. 5. A position velocity diagram of Hi absorption in the cen- 
tral region of the Medusa merger. The absorption has been av- 
eraged over the entire R.A. range of the radio source, where 
the radio continuum is in excess of 0.6 mJy beam -1 . The con- 
tour levels are set at -8, -5.657, -4, -2.282, -2, -1.414, -1 
and 1 times 0.36mJybearrT 1 . The position axis is declination 
relative to 54°48'16'/4. 



by these results, albeit only for a small cross-section through 
the rotating material, but it is apparent that the gas does also 
possess some additional non-circular motions superimposed 
upon this rotation. 

As discussed in the previous section the Hi absorption 
shows significant inhomogeneities in its column densities. In 
particular it is clear that the gas column densities are sig- 
nificantly lower against the bright radio continuum peak in 
the north of the source; along lines of sight where the gas is 
blueshifted relative to the galaxy's systemic velocity. Of course 
these column density variations may, in part, be the result of an 
underlying gas temperature gradient. However such an effect 
would need to be extreme to cause the large observed opacity 
differences and should not considered to be a dominant factor. 
Additionally AH00 find that the distribution of CO emitting gas 
is also asymmetric about the galaxy's kinematic centre. They 
observe bright, compact CO emission components redshifted 
relative to the systemic velocity of the galaxy and offset from 
the radio core by ~ 1"; no equivalent emission features are ob- 
served at similar blueshifted velocities. 

Assuming that the radio continuum emission is related to 
the amount of recent star-formation occurring at anyone point; 
some of the areas with the most active ongoing star-formation 
will be located toward the brightest radio components. Further 
evidence from UV HST observations by Weistrop et al. (2004), 
where high concentration of young (<20Myr old) star clus- 
ters are imaged, confirms that the majority of the recent star- 
formation activity is situated toward the north and western parts 
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of NGC4194. Thus it appears that a considerable portion of 
the activity is situated on the blueshifted side of the source and 
away from the bulk of the cold gas. 

Unlike emission experiments these absorption observations 
only detect gas that is in front of the source. Therefore if the 
star-formation areas to the north and west are being fed by 
the cold gas it is possible that the Hi absorption column den- 
sity toward the radio continuum peaks is depleted. Dynamically 
this might be considered as following the observed rotation ve- 
locity, where the gas approaches the rear of this area of star- 
forming region, is partially converted to stars during its passage 
and hence has a reduced column density in front of this region. 
This scenario implies that the primary reservoir of gas that is 
fuelling these regions of starburst activity is in fact the gas to- 
ward the south of the radio continuum, where both the column 
densities of Hi absorption and CO emission peak (i.e. an area 
along the same line-of-sight as the nuclear dust lane). Of course 
in this hypothesis it is equally plausible that the site of the ac- 
tive star-formation is moving toward the region with high gas 
density as the starburst evolves rather than the gas itself being 
transported into the stationary star-formation region. 

4.4.3. Physics of the cold gas 

As stated in Sec. 3.2 the column density of neutral hydrogen, 
Nh, is related to the gas spin (excitation) temperature, T sp i n and 
the observed Hi absorption by 

N H = 1 .823 x 10 18 T spin j rdV atoms cirT 2 (2) 

where V is the in units of km s - ' and t is the measured opacity 
of the line. Thus any calculated values of Nh will be dependent 
upon the excitation temperature, T sp i n . 

The value of T sp i n for neutral hydrogen is controlled by 
three factors; pumping by Lya photons, the absorption of 
21 cm continuum radiation, and collisions (Field 1959a, 1959b, 
1959c). It is well known that in high density regions that col- 
lisional excitation methods dominate; thus in these regions the 
spin temperature of Hi will approach the gas kinetic tempera- 
ture. For example measurements of Hi emission and absorption 
column densities through dense regions of galactic disks are of- 
ten consistent for excitation temperature in the range of a few 
10s to a few 100 K (e.g. Payne, Salpeter & Terzian 1983; Liszt 
et al. 1983; Braun 1997; Wolfire et al. 2003 and refs therein). 
It should be noted here that these calculations assume that the 
atomic gas, observed by either emission or absorption, is in a 
single component medium in equilibrium and thermalized (i.e. 
nilOOOcm -3 ); see for example Kulkarni & Heiles (1988). The 
assumption that T sp ; n can be approximated to the gas kinetic 
temperature, T&, however, does not necessarily hold true es- 
pecially (if a two component medium is assumed [a cold and 
warm neutral medium]) for areas of lower density gas clouds 
that are subject to irradiation by either significant levels Lya 
photons or 21 cm continuum radiation (Wolfire et al. 2003). 
Work by Gallimore et al. (1999) further investigated these ex- 
citation processes addressing Hi absorbing gas clouds within a 
few tens of parsecs of Seyfert nuclei. One of the primary con- 
clusions of this work was that even in these regions where both 



an abundant 21 cm and Lya source is present, T sp i n will asymp- 
totically approach Tj. for gas densities of nJlOOOctrT 3 (i.e. a 
collision dominated regions). Whereas below this gas density 
the Hi spin temperature, if in the presence of a 21 cm and/or 
Lya photon source, will be significantly influenced by these 
other excitation processes resulting in T sp i n deviating away 
from the kinetic gas temperature (see figure 12 of Gallimore 
et al 1999). However at this point we should also highlight the 
fact that for Hi absorption, t oc (T sp i n AV) _1 . Thus Hi absorption 
studies, such as these, are most sensitive to the coldest atomic 
gas clouds. 

Regarding the Medusa merger we have already assessed 
the plausibility of either an AGN embedded within the central 
starburst of this galaxy or whether the radio emission detected 
can be accounted for by star formation alone (see Sec. 4.3). 
From this discussion, although the presence of a weak or highly 
obscured active core cannot be eliminated, we conclude that 
the dominant energy input into the ISM of the Medusa is via 
the efficient starburst at its centre. Thus this ongoing star for- 
mation will be the primary heating processes for the ISM of 
NGC4194. 

The molecular ISM of the Medusa merger is characterized 
by an elevated CO/ 13 CO 1-0 molecular line ratio of 19 ± 4 
and relatively faint HCN 1-0 emission ( 12 CO/HCN>25) 2 . The 
faintness of HCN emission suggests that the average density 
in the molecular medium is below 10 4 ctrT 3 and we suggest 
that the weak 13 CO 1-0 line is caused by two effects: high 
gas temperatures (Tj >50 K) and the presence of diffuse, un- 
bound gas. High resolution observations reveal that the 13 CO 
1-0 line emission is relatively weak primarily in parts of the 
central dust lane and towards the bright central CO peak (Aalto 
et al in prep). While the latter is likely caused by elevated ki- 
netic temperatures, the former is probably related to diffuse gas 
in a density-wave like situation, similar to that of the bar of 
NGC 7479 (Hiittemeister et al 2000). In general most galaxies 
with elevated values of the CO/ 13 CO 1-0 line ratio also have 
/(60^m//(100//m) Z 0.8 indicating a high average dust tem- 
perature (Aalto et al. 1991, 1995). This further supports the no- 
tion that gas and dust is being heated by the ongoing starburst. 

Considering these arguments it can be concluded that 
the molecular gas within the centre of the Medusa merger 
firstly exists in a relatively low density environment. In this 
sense, although there will clearly be higher density molecu- 
lar gas clouds, the majority of the gas is in regions with n~l- 
5xl0 3 cirT 3 . Additionally the molecular component appears to 
have slightly elevated kinetic temperatures (T& >50K). This 
provides us with an insight into the physics of the molecular 
clouds within the Medusa which can be translated into an up- 
per limit for the density of the atomic gas and a lower limit on 
its kinetic temperature. These limits are valuable in that they 
help to define one extreme of the cold atomic gas phase within 
this merger. The primary caveat is that in most 'cartoon' models 
of the ISM within galaxies the, coldest and densest, molecular 

2 Recent observations by S. Aalto (unpublished) have detected very 
faint HCN emission from the compact, southern molecular peak which 
could be indicative of young or ongoing but obscured star formation 
in this region. 
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clouds are cocooned within an envelope of more diffuse atomic 
gas which will be partially ionised. Thus, apart from the un- 
likely situation where the mixing of the atomic and molecular 
gas is complete, these molecular diagnostic tool can only be 
used to infer limits for the physics of the atomic gas phase. 

5. Conclusions 

Using radio continuum and Hi absorption observations made 
with MERLIN we have studied the neutral hydrogen gas distri- 
bution and kinematics, in addition to the 1 .4 GHz radio contin- 
uum structure of the central kiloparsec of the Medusa merger at 
sub-arcsecond angular resolutions. These observations consti- 
tute the highest angular resolution study, to date, of either the 
neutral or molecular gas within this source. 

Our highest resolution 1 .4 GHz images of NGC 4194 reveal 
a pair of compact mJy radio components situated at the posi- 
tion of the peak in lower resolution radio images of this source. 
These two components are surrounded, primarily toward the 
south by a more diffuse region of weaker radio emission. These 
compact radio components are unresolved with sizes of <25 pc 
and a separation of ~65pc. We have used these radio obser- 
vations to investigate the energy source of the radio emission 
and compare the two most probable scenarios; firstly that some 
portion of the radio emission is related to the presence of a 
weak AGN component in addition to the well known starburst 
emission from this source, or secondly, whether all of the ra- 
dio emission can completely account for the extended nuclear 
starburst. Using a combination of the Hi gas dynamics and star- 
formation rates, derived from the nuclear radio emission com- 
pared, to those derived from observations at other wavelengths, 
we have tentatively concluded that the majority, if not all, of 
the emission is related to star-formation activity. At present we 
cannot categorically eliminate the possibility of a weak AGN 
component without further observational studies such as fur- 
ther high sensitivity multi-frequency VLBI observation or high 
resolution X-ray observations. 

By observing Hi via absorption we have been able to trace 
the gas dynamics in front of the background radio continuum 
within the half central kiloparsec of the Medusa merger. The 
distribution of the absorbing Hi gas shows a significant en- 
hancements toward the south of the radio continuum source. 
This area is along the same line-of-sight as the peaks in the 
12 CO emission observed by AHOO and is also co-spatial with 
a cross-nuclear dust lane (see Fig. 1). Following the results 
of AHOO we conclude that all three of these cold ISM com- 
ponents (dust, CO and atomic hydrogen) are probably phys- 
ically related, with both the CO emission and Hi absorption 
exhibiting similar dynamics. The Hi absorption traces a rela- 
tively shallow, approximately, north-south velocity gradient of 
~320km s _1 kpc 1 over ~ \ kpc of the background radio con- 
tinuum. This gradient is consistent with that observed in CO 
by AHOO. Considering this gradient and the distribution of the 
absorbing Hi and the most active star-forming regions we dis- 
cuss the possible role of the gas dynamics with respect to the 
fuelling of the starburst. Since the Medusa merger has been ex- 
tensively studied via a variety of molecular gas transitions we 
have also discussed the use of these as diagnostic tools to infer 



some limits on the physical environment that the molecular and 
atomic gas resides in. 
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